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The spin configuration induced by single holes and hole pairs doped into stoichiometric, antiferro- 
magnetic cuprates is considered. Unrestricted Hartree-Fock calculations for the three-band Hubbard 
model are employed to study spin-polaron and vortex-like (meron) solutions. Meron solutions for a 
single hole are found to be metastable with higher energy than spin polarons. We observe that the 
meron solution shifts from site-centered to bond-centered as the interaction is increased. Meron- 
antimeron solutions for hole pairs are found to be unstable. The results are in agreement with earlier 
findings for the one-band Ifubbard model. However, we find that the Hubbard interaction of the 
one-band model has to be chosen similar to the one of the three-band model to obtain comparable 
results, not of the order of the charge-transfer gap, as previously expected. 
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I. INTRODUCTION 

Superconductivity appears in cuprates upon introduc- 
tion of electrons or holes into the insulating, antiferro- 
magnetic parent compounds. These are insulators due 
to the strong electron-electron interactions, which lead 
to a Mott-Hubbard splitting of the bands. Upon doping 
with holes, the antiferromagnetic order is lost at about 
2% holes per copper atom and the system becomes super- 
conducting at a somewhat higher doping level. The dop- 
ing of Mott-Hubbard insulators has been studied quite 
extensively. Nevertheless there is still no agreement on 
what happens when a few electrons or holes are intro- 
duced. 

From experiments as well as from band-structure cal- 
culations it is known that the undoped compounds 
have one hole in each Cu d-shell, whereas the oxygen- 
dominated bands are completely filled. Consequently, 
the Cu ions carry spins 1/2, which order antifcrromag- 
netically due to superexchange. Electron and hole dop- 
ing are qualitatively different: Electrons mostly fill the 
Cu hole so that the local spin vanishes. Holes predomi- 
nantly reside in the planar oxygen orbitals, introducing 
another spin 1/2. This physics can be modelled by the 
three-band Hubbard Hamiltonian 

n-iB =edJ2 4a dia + Cp XI P]<y + 

i(T ja 

where (p]<t) is the creation operator of an electron in a 

3d2,2_y2 (2pa;, 2py) orbital at site i (j) and nf^ = dia 
is the corresponding density operator. The parameters 
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ed, Cp denote onsite energies. The hopping amplitude 
between nearest-neighbor copper and oxygen sites (ij) 
changes sign due to the two possibilities of overlap, e.g., 
positive (negative) for the bonds pointing from a Cu site 
in the positive (negative) x and y directions. Finally, Ud 
is the onsite Coulomb interaction for two electrons on the 
same Cu site with opposite spins. 

The Coulomb interaction Up between oxygen electrons 
is irrelevant for the case of one or two holes, since the 
probability of two holes sitting on the same oxygen ion is 
extremely low. Furthermore, the hybridization between 
the oxygen orbitals is omitted as the hopping tpd is the 
dominant process. Furthermore we neglect the Coulomb 
interaction between Cu and O electrons. 

Electron and hole doping also have rather different ex- 
perimental consequences: Antiferromagnetic order is de- 
stroyed by only x 0.02 holes per Cu, whereas it survives 
on the electron-doped side up to a; ~ 0.14. This asymme- 
try is easily understood: Electrons dilute the antiferro- 
magnet, which perturbs the order relatively weakly. On 
the other hand, holes introduce extra oxygen spins, which 
couple antiferromagnetically to both their neighboring 
Cu spins, resulting in a strong ferromagnetic coupling be- 
tween those spins 0. This ferromagnetic bond strongly 
frustrates the antiferromagnet, leading to its rapid de- 
struction Q, 1^ 1^ . This picture is supported by exper- 
iments on cuprates doped with nonmagnetic impurities 
such as zinc nil- They dilute the spin system similarly 
to electron doping, leading to a large critical concentra- 
tion for the destruction of antiferromagnetic order |^, 
which is essentially a percolation transition ^gj . 

The qualitative difference between electron and hole 
doping suggests that a small number of charge carriers 
doped into the parent compounds should also behave 
quite differently. Indeed, several authors propose that 
hole doping leads to the appearance of nontrivial struc- 
tures of the surrounding Cu spins, such as dipolar con- 
figurations P, Q| , vortices or merons 0, I^ J l ^ sk yrm ions 
UM J and domain walls or stripes [§, M, IliL llSl Hal . 
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This should be contrasted with a spin-polaron configu- 
ration, where the spins in the vicinity of the hole might 
be reduced or even inverted in sign but are still coUinear. 
Timm and Bennemann 

PI 

have proposed a mechanism 
for the rapid destruction of antiferromagnetic order in the 
hole-doped regime that relies on the formation of merons 
in the staggered magnetization upon doping. We here use 
the term "meron" to refer to a topological defect in the 
easy-plane Heisenberg model. It is very similar to a vor- 
tex in the XY model, but the order parameter may rotate 
out of the easy plane close to the defect center in a "lo- 
tus flower" configuration. The easy-plane aniso tropy i s 
caused by the Dzyaloshinskii-Moriya interaction ' 
resulting from buckling of the Cu02 planes and spin-orbit 
coupling. For energy reasons, merons are expected to be 
created as meron- antimeron pairs. It is suggested that 
merons appear since each hole introduces a ferromagnetic 
Cu-Cu bond, leading to frustration, which can be reduced 
by forming merons centered on each such bond 0. This 
mechanism only works for hole doping, not for electron 
doping, in agreement with the pronounced asymmetry of 
the phase diagram. 

The question of the existence of meron-antimeron pairs 
in the ground state has been studied with various ap- 
proaches H, 0, 01 . The present situation is not satisfac- 
tory, since unrestricted Hartree-Fock calculations for the 
one-band Hubbard model have found meron-antimeron 
pair solutions, but never as the ground state 0- On 
the other hand, Hartree-Fock calculations have been per- 
formed for the three-band model i20|, but meron- 
antimeron pairs have not been considered. 

The one-band Hubbard model can be derived from 
the three-band model, Eq. (Q, as an effective low-energy 
model |23| . The one-band model is defined by the Hamil- 
tonian 

HiB - E + ii.c.) + [/ ^ cjtcncl^ca. (2) 

In the undoped three-band model the Fermi energy lies 
between an occupied oxygen-type band and the unoccu- 
pied upper Hubbard band. The charge-transfer gap be- 
tween them is mimicked by the Mott-Hubbard gap of the 
one-band model and the onsite interaction U of the one- 
band model must be chosen accordingly |2^ |2^ . How- 
ever, this derivation 21] of the one-band model assumes 
a polaron-like configuration for a doped hole and would 
not work if meron-antimeron pairs were the true ground 
state. 

On the other hand, a slave-boson approach did find 
meron-antimeron pairs as the ground state for two holes, 
again in the one-band model '9] . This approach goes be- 
yond the Hartree-Fock approximation (HFA) and leads 
to lower variational ground-state energies. Furthermore, 
in a one-band spin-flux model Berciu and John also 
find meron-antimeron pairs as the ground state, employ- 
ing the configuration-interaction (CI) method, which also 
goes beyond Hartree-Fock theory. Both approaches have 



not yet been applied to the three-band model. The ques- 
tion appears to be open. 

In the present paper we consider the three-band Hub- 
bard model within the HFA. In Sec. we briefly out- 
line the theory. In Sec. IIIII we compare the results with 
Hartree-Fock calculations for the one-band model for 
larger system sizes than studied previously. We study the 
stability of spin-plasmon and meron solutions for a sin- 
gle hole and, in particular, of meron-antimeron solutions 
for two holes. We also consider changes in the meron 
structure as a function of the interaction strength. 



II. THEORY 

Since it is mostly standard, we can be brief in describ- 
ing our implementation of the unrestricted Hartree-Fock 
approximation 24j. We here do this for the three-band 
Hubbard model. The HFA reduces the many-body prob- 
lem to a set of single-electron problems by decoupling the 
interaction, in our case in the particle-hole channel. To 
include transverse spin degrees of freedom we decouple 
on the level of creation and annihilation operators, 

(3) 

with the spin-flip operators S""*" = + iS^ = and 
= S"^ — iS"^ = (i|c?| and the constant energy shift 

Eo = U,J2 (li^f)' - (Sn' - {Sff - {Sf)"] . (4) 

i ^ ' 

After decoupling the interaction, the many-particle 
ground state \^) for T = in the HFA is constructed by 
creating electrons in the lowest-energy single-particle 
states out of the vacuum |0), 

We 

i*)-n4io>. (5) 

The single-particle states are given in terms of wave func- 
tions (/)„ by 

(^\^^4>-a{i,(y)c\„, (6) 

where the sum over i runs over all relevant Cu ii^i^yi and 
oxygen p^; and Py orbitals, and cj^ creates an electron in 
orbital i. 

The Hartree-Fock eigenequations are derived by min- 
imizing the expectation value to the total energy in the 
state l^*) with the cr) as variational parameters. For 
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the n-th state with spin cr =|, J, at the Cu site i the equa- 
tion is of the form 

En Ipnih cr) = -t^d X! ^) + 

The notation j G means that the sum is performed 
over the sites j which are nearest neighbors of site i. 
cr = ((T^, (7^, cr^) is the vector of Pauh matrices and 
is the unit matrix. The local Hartree-Fock fields of the 
charge and spin density, 

{<) = mY.dlAM'^) = Y.\^n{i,^)?, (8) 

fiiy 

= ^'YKihtJ')f^tJ^i^'^>nii,'^) (9) 

are computed self-consistently. The sum over n is per- 
formed for the A^e electrons in the system. The terms for 
the oxygen sites j are of an analogous form. Finally, the 
Hartree-Fock ground-state energy is given by the sum of 
the eigenvalues £"„ for occupied states, 

n=l 

We solve the Hartree-Fock equations JT)) by iteration, 
which requires a numerical diagonalization at each step. 
Self-consistency is assumed when the difference in all spin 
and charge components is less than e ~ 10~^. Typically, 
the charge distribution converges very rapidly to the self- 
consistent value, while the spin configuration needs a 
lot more iteration steps. The number of steps depends 
strongly on the initial configuration and the boundary 
conditions. 



III. RESULTS AND DISCUSSION 

In the present section we examine the Hartree-Fock so- 
lutions for the three-band Hubbard model for one and 
two holes (2J|. The results are compared with corre- 
sponding ones for the one-band model. The calculations 
have been carried out for a square lattice with up to 
14 X 14 unit cells for the three-band model as well as for 
the one-band model. We have treated periodic as well as 
open boundary conditions (BC). The latter disturb the 
elementary cells at the edge of the lattice in the three- 
band case. We have studied edges of alternating copper 
and oxygen ions but also edges made up entirely of oxy- 
gen ions. Finite size effects are stronger for open BC, 
as expected, but are neglectible for the lattice sizes used 
here. Relevant parameters in the three-band model are 



(i) the hopping tpd, (ii) the Coulomb interaction Ud for 
copper sites and (iii) the band gap A = — between 
the oxygen and copper orbitals. A Coulomb interaction 
between oxygen electrons and a non-zero oxygen-oxygen 
hopping amplitude do not change the solutions qualita- 
tively and are neglected in the present study. A band gap 
A > in combination with C/d > A implies that holes 
preferentially occupy oxygen orbitals as observed in ex- 
periments. In this regime the occupied oxygen states 
and the empty copper states in the upper Hubbard band 
are separated by a charge-transfer gap, which is given by 
[/(J — A. For the calculations we have used the parameters 
taken from Ref. 

ipd = leV, C/rf = 8eV, A = 5eV. (11) 

The results for the three-band model are compared 
to corresponding Hartree-Fock results for the one-band 
Hubbard model with t = 1 eV and U in the range 
from 5 to 8 eV. 
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FIG. 1: Self-consistent Hartree-Fock solution of the three- 
band Hubbard model for one hole. The figure shows part 
of a 10 X 10 lattice with the spin and relative charge den- 
sity indicated for each atomic site by the arrows and circles, 
respectively. The reduced copper spin at point (5, 5) forms 
a spin polaron with its four nearest copper neighbors. The 
charge is mainly distributed over the four oxygen neighbors. 



A. Single-hole solutions 

For a single hole we find a spin-polaron configuration 
as the self-consistent solution of the HFA for open as 
well as periodic boundary conditions, see Fig. ^ In its 
center a reduced copper spin (of modulus 0.55 in units 
where the fully polarized state corresponds to unity vs. 
an average order parameter of 0.78) is pointing in the 
same direction as its four nearest copper neighbors, e.g.^ 
a ferromagnetic "microdomain" is formed. The charge 
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is mainly distributed over the four oxygen neighbors of 
the central copper ion. Small magnetic moments at the 
oxygen sites are induced, which are oriented oppositely to 
the copper spins. This result is consistent with the over- 
compensation of the antiferromagnetic superexchange by 
an effective ferromagnetic Cu-Cu coupling for all values 
of Ud > A, supporting the picture of Aharony et al. 0. 
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FIG. 2: One-particle energy spectra for self-consistent solu- 
tions of the HFA. LHB and NB denote the occupied lower 
Hubbard and non-bonding bands, respectively. UHB denotes 
the unoccupied upper Hubbard band |23. |23(I . (a) Spectrum 
for the spin-polaron solution of Fig. (b) Spectrum for the 
meron solution of Fig.|3Ia). (c) Spectrum for the solution of 
Fig. |S| The right-hand-side figures are enlargements of the 
regions close to the Fermi energy. 

In the one-particle spectrum of Fig. |2Ia) two states 
appear in the gap due to the interaction. One is split 
off from the occupied oxygen-dominated band but is now 
empty due to the extra hole. This is the Zhang-Rice 
singlet lij. The other state originates from the unoccu- 
pied upper Hubbard band. The first unoccupied state is 
to about 75% localized on oxygen sites emphasizing the 
oxygen character of doped holes '2^ . 

It is worth considering the limit of str ong interaction 
Ud/tpd > 1 studied by Zhang and Rice [Ij. For Ud — 
50 eV, A = 25 eV wc find a polaron solution similar to 
Fig. ^ but with the copper spin moment in the polaron 
center hardly reduced (to 0.99) and the charge mainly 
(to 80%) localized on its four oxygen neighbors. These 
four ions together carry a spin moment of —0.79. Thus 
in the HFA the central copper spin forms a singlet with 
its four oxygen neighbors, as predicted . 

It is instructive to compare the results with the one- 
band Hubbard model, for which the polaron configura- 
tion is also a self-consistent solution in the HFA 8]. For 



an intermediate interaction of t/ ~ 5 ... 8 eV the charge 
distribution and the spin moment at the polaron center 
are of the same order of magnitude as in the three-band 
case so that in this respect the one-band model indeed 
mimicks the three-band physics. However, it is surprising 
that we find this relation for values of U for the one-band 
model similar to Ud = 8eV and not to the charge-transfer 
gap Ud — A = 3 eV, as we would have expected j2.2, 23]. 
On the contratry, for U < 3eV we find the diagonal cigar- 
shaped configuration predicted in Refs. [23, without 
any inverted Cu spin We suggest that U is of the 

order of Ud since the local moment is dominated by the 
interaction in the copper d orbital and is only marginally 
affected by the hybridization with oxygen orbitals. 

We now turn to meron configurations for a single doped 
hole. We start the HFA with a coplanar vortex con- 
figuration. Site-centered Merons turn out to be stable 
Hartree-Fock solutions for Ud 7 eV. Because of the 
non- vanishing winding number only open boundary con- 
ditions are possible. In Fig. Ofa) a site-centered meron 
solution for Ud = 7 eV is shown for a lattice with edges 
of alternating copper and oxygen ions. But there is no 
qualitative difference to lattices with edges made up en- 
tirely of oxygen ions. Increasing Ud shifts the center of 
the meron first in the direction of a plaquette and then 
onto a bond, see Figs.OJb) and (c). In addition, a ten- 
dency towards ferromagnetic alignment, i.e., formation 
of a spin polaron, is observed. These results correspond 
to the one-band solutions for [/ > 8 eV, in which the 
spin in the meron center is no longer suppressed and 
shows a tendency towards formation of a ferromagnetic 
microdomain. Note that again U has to be chosen sim- 
ilar to Ud and not to the charge-transfer gap to obtain 
this result. To exhibit the transition from site- to bond- 
centered merons more quantitatively, we plot in Fig. ^ 
the average (S)^ at the site i = (6, 6) for the three-band 
solution and the corresponding quantity for the one-band 
model. For a site-centered meron this number is negligi- 
ble due to symmetry. 

For the one-band model we have compared the ener- 
gies of polaron and meron solutions for a single hole by 
considering their respective formation energies As 
expected, the meron energy is higher than the energy of 
the spin polaron and diverges logarithmically with sys- 
tem size. However, the iteration process of the HFA 
turns out to be unable to change the winding number. 
Adding small random numbers during the iteration pro- 
cess makes it converge to a spin-polaron configuration. 
Thus a meron is a metastable solution. 

In the one-particle spectrum for the meron solution 
with Ud ~ 7 eV, Fig. [^Jb), two discrete, degenerate hole 
states appear: one from the filled oxygen-dominated non- 
bonding band and the other from the unoccupied upper 
Hubbard-band. The states are degenerate because an 
additional electron can be introduced either with spin up 
or spin down. The same holds for the one-band model. 
A mirror image of these degenerate states appears in the 
gap between the occupied bands. 
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FIG. 3: Selfconsistent Hartree-Fock solutions with (a) Ud — 
7eV, (b) Ud = 8eV and (c) Ud = lOeV on an 11 x 11 lattice. 
For Ud < 7eV site-centered merons are stable solutions of the 
HFA. Increasing Ud shifts the center of meron but preserve 
the winding number of the initial configuration. 



B. Two- holes solutions 

We now turn to a system doped with two holes to in- 
vestigate possible meron- antimeron solutions. As initial 
configurations we consider pairs of bond-centered as well 
as plaquette-centered defects, see Fig. El and both open 
and periodic BC. We find that these configurations are 
not stable under the HFA iteration. This result holds 
for all parameter sets we have investigated. Instead, the 
system develops into a polaron-type spin configuration, 
the detailed structure of which depends on the choice 
of BC. For example. Fig. shows the solution for the 
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4 6 8 10 12 

Ud, U (eV) 

FIG. 4: The average value (S)^ at the site closest to the meron 
center as a function of the interaction U for the one-band 
model and Ud for the three-band model. For a site-centered 
meron, (S)^ is negligible due to symmetry. A non-vanishing 
value corresponds to a qualitative change in the spin texture, 
see FigEl 



bond-centered initial configuration and periodic BC on a 
14 X 14 lattice. We find two spin polarons with antifer- 
romagnetically aligned center spins. The polaron centers 
are localized on neighboring copper sites. At the oxy- 
gen sites small magnetic moments are induced which are 
oriented in antiparallel to the copper spins. 

The origin of the instability of meron-antimeron solu- 
tions is their higher energy and their vanishing winding 
number. The latter means that they can be destroyed 
during the HFA iteration since no high energy barrier has 
to be overcome. Since for large separations R of meron 
and antimeron the interaction increases logarithmically 
with R, the energy of the meron-antimeron solution will 
eventually become larger than the spin-polaron energy. 
However, our results indicate that this is also the case at 
the smallest separations. 

Again, similar results are obtained for the one-band 
model. Bond-centered as well as plaquette-centered 
meron-antimeron pairs as initial configurations develop 
into spin polaron solutions as shown earlier by Verges et 
al. 0. In contrast, Scibold ^] has found a stable solution 
in the one-band model using a slave-boson approximation 
beyond Hartree-Fock. 

In Fig. |2Ic) the spectrum of the discussed double- 
polaron solution is shown. For each doped hole one state 
splits off from the occupied non-bonding band and an- 
other from the unoccupied upper Hubbard band. The 
first two hole states are localized to about 80% on oxy- 
gen sites. Further, the states are energetically degenerate 
reflecting the symmetry of the solution under inversion in 
the center of the bond between (7, 6) and (7, 7) in Fig. El 
and simultaneous spin inversion. 
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FIG. 5: Starting configurations for tfie spin and cliarge den- 
sity for a system with two holes on a 14 x 14 lattice: (top) 
plaquette-centered meron-antimeron pair and (bottom) bond- 
centered meron-antimeron pair. 

IV. CONCLUSIONS 

In this paper we have used unconstrained Hartree-Fock 
theory at T = 0, keeping all local spin and charge de- 
grees of freedom to study inhomogenous spin textures in 
the three-band and the one-band Hubbard model. This 
approach neglects flucutations around the mean-field so- 
lutions. Thus one has to consider the reliability of the 
these solutions. The HFA predicts the correct behavior 
for large Ud in the half-filled case. For one doped hole 
spin polarons are obtained as solutions for the three-band 



Hubbard model. This agrees with the case of the one- 
band model (Si] , in agreement with the mapping |^ 
of the three-band on the one-band model at low energies. 
However, we find that in the one-band model an on-site 
interaction U of the order of the three-band interaction 
\Jd is required to obtain similar results. 
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FIG. 6: Resulting double-polaron configuration for two holes 
on a 14 X 14 lattice with periodic boundary conditions using 
the parameters in Eq. IIIH . 



A meron is only a metastable self-consistent solution 
due to its nonzero winding number. As our main result 
we find that meron-antimeron pairs are not stable solu- 
tions in the HFA of the three-band Hubbard model. This 
holds for periodic as well as open boundary conditions. 
The same is true for the one-band model. Thus a consis- 
tent picture of the one- and three-band models emerges, 
at least in the framework of the HFA: Both show the 
same qualitative behavior for weak hole doping. 

While the mean-field approximation serves as an use- 
ful tool for studying the local charge and spin density, it 
yields only an upper bound for the true ground-state en- 
ergy. In the one-band Hubbard model Seibold Q found 
a meron-antimeron ground state for T = for suitable 
parameters using a slave-boson method. This method 
uses a Hartree-Fock ansatz as the initial configuration 
but treats the spin and charge excitations at each site as 
distinct bosonic degrees of freedom. It leads to a lower 
ground state energy than the HFA and is thus the su- 
perior variational method. We suggest that this method 
applied to the repulsive three-band Hubbard model could 
clarify the still open question of stable meron-antimeron 
pairs as the effect of hole doping. 
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